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Abstract. Strong broad emission lines are the most important signatures of active galactic 
nuclei. These lines allowed to discover the cosmological nature of quasars, and at present these 
lines allow for convenient method of weighting the black holes residing in their nuclei. However, 
a question remains why such strong lines form there in the first place. Specifically, in the case of 
Low Ionization Lines, there must be a mechanism which leads to an efficient rise of the material 
from the surface of the accretion disk surrounding a black hole but at the same time should 
not give a strong signature of the systematic outflow, as the Balmer lines are not significantly 
shifted with respect to the Narrow Line Region. We determine the effective temperature of the 
accretion disk underlying the H/3 line at the basis of the time delay measured from reverberation 
and the simple Shakura-Sunyaev theory of accretion disks. We obtain that this temperature 
is universal, and equal 995 ± 74 K, independently from the black hole mass and accretion rate 
of the source. This result suggests to us that the dust formation in the disk atmosphere is 
responsible for the strong rise of the material. However, as the material gains height above 
the disk it becomes irradiated, the dust evaporates, the radiation pressure force suddenly drops 
and the material fall back again at the disk. Therefore, a failed wind forms. In the simple 
version of the model the disk irradiation is neglected, but in the present paper we also discuss 
this irradiation and we use the observed variation of the Broad Line Region in NGC 5548 to 
constrain the character of this non-local non- stationary phenomenon. The current instruments 
cannot resolve the Broad Line Region but future instrumentation may allow to test the model 
directly. 



1. Introduction 

Broad emission lines in quasars and in other types of active galaxies have lead to the identification 
of those objects as distant sources with strong non-stellar nuclear emission coming from their 
nuclei due to accretion onto black holes. Both Carl Seyfert (1943) pQ and Maarten Schmidt 
(1963) [2] in their seminal papers were just studying emission lines. Since that time there was a 
tremendous progress in our understanding of the nature of active galaxies and in our knowledge 
of the Broad Line Region (BLR) properties. Emission line measurements are also nowadays a 
powerful tool to measure the mass of the central black hole (Kaspi et al. [3j; Peterson et al. [4]) 
and these measurements done for large samples of quasars (e.g. Vestergaard et al. [5]) give the 
insight into the evolution of galaxies. But one question still remains to be answered: why this 
region actually forms there. 

The properties of the BLR contain a hidden hint. The argument goes through three 
independent lines of reasoning. 



First point: The lines cover a considerable fraction (0.1 - 0.3) of an AGN sky from the point 
of view of the nucleus, as estimated from the line luminosity (ref. [6]). At the same time they are 
almost never along the line of sight to the nucleus although NAL are seen in a significant fraction 
of objects (BAL are rare; 20% of normal quasars [Tj) and both NAL and BAL form at larger 
distances from the nucleus. There are rare reports of the transient absorption phenomenon seen 
in X-ray data which are likely to be due to BLR clouds [8] . These observations indicate that the 
BLR is not spherical but flat and located close to the equatorial plane; in side-viewed objects 
this region is hidden from us by the dusty-molecular torus and in top- viewed objects the line of 
sights does not cross BLR. Kinetic arguments also suggest such a flattening (e.g. Nikolajuk [S], 
Collin IIOJ). 

Second point: the line widths are from 1000 to 10 000 km s~^, consistent with the Keplerian 
disk-like motion. However, lines rarely show two-component profile expected from a flattened 
disk-like configuration. There are rare objects showing clearly double-peaked lines but these 
are mostly radio-galaxies (Eracleous & Halpern |11]), and some objects hint for the disk-like 
profiles in the variable part of their line profiles. The majority of quasars have a single peak. 
Murray & Chiang have shown that even a disk-like configuration can lead to single profiles but 
it requires special wind velocity profile [T2]. And no systematic shift (outflow velocity signature) 
with respect to the Narrow Line Region (NLR) is seen in Balmer lines. 

Third point: The careful studies of the luminosity and the line ratios, particularly in respects 
to the Low Ionization Lines like Balmer lines and Mg II (Collin-Souffrin [l3]), suggest that 
additional collisional heating is needed in addition to pure radiative heating estimated from the 
shape of the continuum emission from the central region of an active galaxy. 

All three problems point towards the BLR as flat configuration roughly in Keplerian motion 
but with additional turbulence with the turbulence velocities of order of a thousand km s~^. 
Such velocities, although only a fraction of the Keplerian speed would be 100 times higher 
than the sound speed, and the development of a such highly supersonic turbulence is not easily 
explained. 

In our paper we offer the plausible mechanism which gives strict quantitative predictions 
where the LIL part of the BLR would form. Future high resolution observations will possibly 
resolve the BLR region and offer direct test of the proposed mechanism. 

2. Disk temperature underlying the BLR 

We determine the effective temperature of the accretion disk atmosphere underlying the BLR 
using the observational determination of the BLR radius from the reverberation method and 
combine it with the Shakura-Sunyaev theory of accretion disks |14j . 

The scaling of the BLR with luminosity discovered by Kaspi et al. [3] has been studied by 
several authors (e.g. Peterson et al. [3]). Here we use the result of Bentz et al. [15], where the 
starlight contamination was carefully removed. The power-law index in the scaling between the 
BLR size and the monochromatic luminosity in their analysis is consistent with 0.5 within the 
errors. We thus fixed this slope at 0.5, and refitted their data again, obtaining the relation 

logRBLRim = 1-538 ± 0.027 + 0.5 log ^44,5100, (1) 

where -Rbl/?[H/3] is in light days and ^44,5100 is the monochromatic luminosity at 5100 
A measured in units of 10^'^ erg s^^. The value 0.027 is the error in the best-fit vertical 
normalization; the dispersion around the best fit is larger but still small at 0.21. 

We assume that the disk is not strongly irradiated, and the outflow from the disk is not too 
strong, so the simple theory of the alpha-disk (Shakura & Sunyaev [13]) applies. In this case, 
the monochromatic flux at 5100 A can be calculated from the model, if the black hole mass 
and the accretion rate are known. Taking the numerical formula of Tripp et al. [16] (corrected 



for a factor-of-two error, Nikolajuk, private communication), we derive 



log L44,5ioo = ^ log(MM) - 43.8820 + log cos i, (2) 

where the black hole mass, M, and the accretion rate, M, are in g and g s~^, respectively. This 
formula allows us to calculate the BLR radius from the Eq. [H if the product of (MM) is known. 

For completeness, we introduce a geometrical factor representing the effect of the inclination 
angle, i for the measured delay (for details, see Czerny &: Hryniewicz [17j). This gives the 
relation between the actual disk radius, r, and the measured delay of H/S line: 

, = RBLRim 
1 + sin i 



and we fix this factor adopting the average inclination i = 39.2° after Lawrence & Elvis 

The Shakura-Sunyaev accretion disk theory allows us to calculate the effective temperature 
of the disk at that radius as a function of the black hole mass and accretion rate 



where ctb is the Stefan-Boltzmann constant. Since BLR is far from the central black hole, the 
inner boundary effect can be neglected. 

When we combine Eqs. [H O El and HI the dependence on the unknown mass and accretion 
rate vanishes, and we obtain a single value for all the sources in the sample 

Te/j = 995 ±74K, (5) 

where the error reflects the error in the constant in Eq. [TJ We stress that the resulting value 
does not depend either on the black hole mass or accretion rate. The value is universal, for all 
sources. 

This value is interesting since it is close to the critical temperature at which the dust can 

1/2 

form. The distance is smaller than the distance to the dusty torus Rdust ~ OAL^^ pc (Nenkova 
et al. |19j : where -L45 stands for the bolometric luminosity in 10^^ erg s~^) since in our case the 
temperature results only from the local disk dissipation instead of the irradiation by the central 
parts of AGN. 



3. The mechanism of formation of the BLR 

This uniform effective temperature of 1000 K of the disk underlying the onset of the BLR implies 
that the dust present in the disk atmosphere is responsible for the phenomenon. The radiation 
pressure acting on dusty atmosphere is very large, as it is known for example for Wolf-Rayet stars 
and the strong wind develops. However, unlike in stars, the further development of the wind 
is modified by two factors: gravitational attraction in disk geometry rises with the height and, 
more importantly, the irradiation by the central parts of the accretion disk (UV) also rises very 
strongly due to better geometrical exposure. The dust evaporates, the rising material looses the 
radiation pressure support and falls back onto the disk surface. Such a failed wind can account 
both for the large additional turbulent velocity without any strong signature of outflow since 
both outflow and inflow take place with comparable velocities, and for the additional heating, 
as the outflowing and infalling clumps collide. 

This region is separate from the region known as dusty /molecular torus. At large distances 
the UV irradiation is less important, the dust survives even at large height from the disk and 
the outflow proceeds forming an appearance of a torus. The schematic representation of the 
BLR formation mechanism is shown in Fig[TJ 
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Figure 1. The schematic picture of the 
formation of BLR. In this region dusty 
wind rises up, the dust evaporates due to 
irradiation by the central parts, gravity wins 
and dustless material falls back onto the 
disk. Dusty torus is further from the nucleus 
where irradiation does not destroy the dust. 



Figure 2. The dependence of the size of the 
BLR in NGC 5548 on the monochromatic 
luminosity: data points come from Bentz 
et al. [15], and the best fit line has 
the slope 0.90 it 0.22 in the logarithmic 
space, instead of 0.5 characteristic for the 
statistical sample. 



4. Discussion 

The mechanism of BLR formation proposed in this paper and in Czerny &: Hryniewicz [17J is 
attractive since it is universal for all sources, and the physics of dusty winds is known from 
stellar evolution. Although the possible presence of the dust in the BLR clouds was previously 
discussed by a number of authors, our idea of the key importance of the dust for BLR is new 
and should be further tested. The whole idea automatically explains why the size of the BLR 
region scales with the square root of the monochromatic luminosity, with small dispersion, since 
no dependence on the second parameter - for example bolometric luminosity is expected. We 
expect a weak dependence of the inclination angle but this also remains to be tested. 

4.1. Determination of the disk temperature in all monitored sources 

In the paper by Czerny & Hryniewicz [17j we tested additionally the mechanism by determining 
the disk BLR temperature for a number of sources which were carefully monitored for many 
years (e.g. Peterson et al. ^j) and their luminosity was corrected for the starlight contribution 
at the basis of the HST data (Bentz et al. [15j). For all sources we obtain the temperature value 
close to 1000 K. This is also the case for NGC 5548 monitored for 14 years. 

We do not see any correlation between the accretion rate and the disk effective temperature 
at the onset of BLR in the objects for which the reverberation monitoring was performed. 
Mechanism simply connecting the BLR distance with irradiation would likely contain the trend 
with the bolometric luminosity which (in Shakura-Sunyaev model) rises faster with the accretion 
rate than the monochromatic luminosity. 

4-2. The role of the disk irradiation 

In our approach to the description of the disk atmosphere underlying the BLR and to the 
description of the monochromatic flux we used the simple local stationary disk model of Shakura 
& Sunyaev. However, observations of AGN show that AGN are variable and this variability is 



a central point of the monitoring approach and the measurement of delays. It means not only 
that the disk is not stationary, but it implies that non-local phenomena (irradiation of the outer 
disk parts by the inner disk) are complicated by these time delays. 

The measured time delays between the continua measured at UBV bands are typically of 
order of a day or less in Seyfert galaxies (see e.g. Sergeev et al. [20]), the major trends in optical 
luminosity are of order of tend-hundreds of days, and measured time delays of H/3 line are of 
order of days (see compilation in Bentz et al. P[S]). Both an outer and an inner disk are the 
sources of variability, as argued for example by Gaskell |21j . and the perturbations in the outer 
disk propagate mechanically inwards through the disk (sound waves, viscous propagation) while 
perturbations in the inner disk affects the outer disk through irradiation, and the timescale is 
the light travel time. This complex situation can be roughly represented with two different 
accretion rates/luminosities: inner accretion rate visible through the monochromatic flux and 
responsible for the current strength of irradiation and the outer accretion rate (past value, likely 
closer to the average value) characteristic for the disk structure underlying the BLR. 

We can now consider various possibilities depending on the importance of the irradiation on 
the measured value. 

If the whole idea of Czerny & Hryniewicz [T7j is incorrect, and the dust formation in the 
accretion disk atmosphere has nothing to do with the formation of the BLR, then we return 
to the classical point of view that the onset of BLR requires a specific value of the ionization 
parameter. In this case the radius of the BLR should scale with the accretion rate at the inner 
disk as 

RbLR OC Ml/2. (g) 

If we assume that the monochromatic flux is nevertheless dominated by the local dissipation, 
i.e. we can use Eq. [21 the final dependence between the monochromatic luminosity and the BLR 
radius will be 

Rlbr cx LOfooM"'/', (7) 

i.e. the slope will be close to the older value and the presence of the additional dependence 
on the black hole mass would be expected. For a given source, only the dependence on the 
monochromatic luminosity matters, while in a sample of sources there is a likely bias towards 
positive correlation between the monochromatic luminosity and the mass, leading to much more 
shallow dependence. 

Therefore, we additionally plotted the results of the monitoring for NGC 5548, taking the 
data points from Bentz et al. where the starlight was subtracted. The relation is shown in Fig.[2j 
The slope is indeed steeper, equal to 0.90 ± 0.22. The lower limit is consistent within the error 
with the expectations based on the ionization parameter, but the value itself is surprisingly high. 
Therefore, we considered whether the dust based scenario can explain such a steep slope. If we 
follow the general idea of Czerny & Hryniewicz [17j of dust onset of the BLR but we include 
the irradiation of the disk surface we also obtain the dependence of the BLR position on the 
accretion rate at the inner disk. Assuming that the irradiation actually dominates we have 

F,.s, = ^-^fi^l (8) 

where h{r) is the disk thickness, and the radial dependence is in the form of a power law, 
h{r) OC r'^. 

Loska et al. |22] considered direct and indirect irradiation of an AGN accretion disk and 
calculated the disk thickness of the non-irradiated and irradiated disk taking into account the 
vertical disk structure with convection and realistic opacities, with atomic transitions, molecules 
and dust. The disk thickness rises slowly with the radius at the inner part of the disk and then 
flares at larger distances, roughly as predicted by the original Shakura & Sunyaev model. The 



Table 1. Observed size of the BLR in the monitored AGN. 
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final rise in Loska et al. [22] is roughly with h oc r^ "^^ (see their Fig. 3). If we supplement this 
value of 7 into Eq. [HI assume the BLR location at a fixed temperature determined by the dust, 
and combine it with the expression for the monochromatic fiux we obtain the expected scaling 

Rblr oc Lgfo>^"°-'^ (9) 

i.e. the trend with the black hole mass appears but the trend with the luminosity matches nicely 
the one seen for NGC 5548. 

Instead, if the irradiation dominates at 5100 A the spectrum has much flatter slope 
(-Fjy cx than in the Shakura-Sunyaev disk {Fy oc z^^^'^), and the monochromatic flux is 

proportional to the accretion rate, i^sioo oc M. In this case the combination of the irradiation 
dominating 5100 A and the assumption of the speciflc value of the ionization parameter specified 
BLR gives the required relation for the statistical sample 

Rblr cx Lgfoo> (10) 
but it is inconsistent with the trends observed in NGC 5548. 

^.3. Further tests of the model 

Better tests of the idea that the dust is the key mechanism behind the formation of the BLR 
require both more results of the reverberation monitoring and better description of the actual 
dusty outflow. 

Observationally, results from many years for a single source would allow to see whether the 
steeper relation seen for NGC 5548 is a typical pattern for a single source variability. Monitoring 



of more objects in a still broader range of masses and accretion rates would allow to test the 
exact form of coupling of the BLR size to the mass and accretion rate. 

Direct observational approach to the location of the BLR is not possible at present. The 
requested spacial resolution is very high: the optimum source is NGC 4051 (100 /^arcsec) due to 
its low redshift (see Table [1]), but in the future even such an extreme resolution will be achieved. 
Such a possibility was already envisioned by Elvis & Karovska |23j . 

Theoretical studies are also needed since the actual formation of the dust in the disk 
atmosphere, dusty outflow, dust/gas coupling, dust evaporation and subsequent inflow are very 
complex phenomena. There are papers addressing some of the issues (for example, see a recent 
paper of Dorodnitsyn et al. [23] about the dusty outflow forming the dusty torus), but never all 
of them. In addition, the dust may acquire charge, and the presence of the large scale magnetic 
field may also modify the dust particle behaviour, (see e.g. Kovar et al. for extreme examples 
of charges close to a black hole). 
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